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Abstract 

We analyze the decays {h°, H°, A } — ► sb within two Higgs Doublet Models with 
Natural Flavor Conservation (2HDM) type I and II. It is found that the Higgs bosons 
decay into bottom-strange can lead to a branching ratio in the range 10 -5 — ► 10~ 3 for 
small tan/3 ~ 0.1 — * 0.5 and rather light charged Higgs in the 2HDM type I. When 
tan/? > 1, one can easily reach a branching ratio of the order 10 -5 . In 2HDM type 
II, without imposing b — > S7 constraint, the situation is the same as in 2HDM type I. 
If b — > 57 constraint on charged Higgs mass {Mh± > 350 GeV) is imposed, we obtain 
Br(h° — > sb) in the range 10- 5 -10" 6 . A comparison between the rates of h° — > sb 
and h° — > 77 is made. It is found that in the fermiophobic scenario, h — > 77 is still 
the dominant decay mode. 
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1 Introduction 



One of the goals of the next generation of high energy colliders, such as the large hadron 
collider LHC [1] or the linear collider LC [2] or muon colliders, is to probe top Flavor- 
Changing Neutral Couplings 'top FCNC as well as the Higgs Flavor- Changing Neutral 
Couplings 'Higgs FCNC. FCNC of heavy quarks have been intensively studied both from 
the theoretical and experimental point of view. Such processes are being well established in 
the Standard Model (SM) and are excellent probes for the presence of new physics effects 
such as Supersymmetry, extended Higgs sector and extra fermions families. 

Within the SM, with one Higgs doublet, the FCNC Ztc vanishes at tree-level by the 
GIM mechanism, while the ytc and gtc couplings are zero as a consequence of the unbroken 
SU(3) C x U(l) em gauge symmetry. The Higgs FCNC Hie and Hsb couplings also vanish 
due to the existence of only one Higgs doublet. Both top FCNC and Higgs FCNC are 
generated at one loop level by charged current exchange, but they are very suppressed by 
the GIM mechanism. The calculation of the branching ratios for top decays yields the SM 
predictions [3], [4]: 

Br(t -> Zc) = 1.3 x 1CT 13 , Br(t -> 7 c) = 4.3 x KT 13 , Br(t -> gc) = 3.8 x 1CT 11 , 
Br(t -> He) = 5.6 -> 3.2 x KT 14 for M H = 115 -> 130 GeV. (1) 

While for Higgs FCNC, calculation within SM leads to: 

Br(H -> sb) « 1(T 7 (resp 10~ 9 ) m H = 100 (resp 200) GeV 

Br(H -> tc) « 1.5 x 10~ 16 (resp 3 x 10~ 13 ) m H = 200 (resp 500) GeV (2) 

Many SM extensions predict that these top and Higgs FCNC can be orders of magnitude 
larger than their SM values (see [5] for an overview). For the Higgs FCNC, an important 
class of models where Higgs FCNC appear at tree level are the so called Two Higgs Doublet 
Model without Natural Flavor Conservation (NFC) 2HDM-III [6, 7, 8, 9]. In this class of 
models, the branching ratio of h — > tc can be larger than 10% in some cases [7]. In the 
framework of 2HDM with NFC type I and II, top and Higgs FCNC have been studied in 
[10, 11]. It was shown that in 2HDM-II the £r($ -> tc), $ = h° or H°, may reach 10" 5 
for CP-even states [11]. This rate is almost eight orders of magnitude larger than the SM 
one. 

Top and Higgs FCNC couplings have been addressed also in supersymmetry [12, 13, 14, 15, 
16]. In those studies it has been shown that Br(h° — > sb) can be in the range of 10 _4 -10 -3 . 
This rate originates mainly from the flavor violation interactions mediated by the gluino 
[12, 14]. In case of MSSM with R parity conservation, the top FCNC coupling t — > ch°, 
can reach 10~ 4 branching ratio [15] in case of flavour violation induced by gluino. 

Hence, Higgs and top FCNC offer a good place to search for new physics, which may 
manifest itself if those couplings are observed in future experiments such as LHC or LC 
[1, 2]. Therefore, models which can enhance those FCNC couplings are welcome. 

The aim of this paper is to study Higgs FCNC couplings such as $ — > sb, $ = h°,H°,A°, 
in the framework of NFC two Higgs Doublet Models type I and II. It is found that the 
branching ratios of Br(& — > sb), $ = h°,H°,A°, can be greater than > 10~ 5 in quite a 



2 



substantial region of the 2HDM parameters space. Br($ — > tc) requires large tan/5 and 
light charged Higgs [11] while £?r($ — > s6) requires rather small tan/5 together with light 
charged Higgs and large soft breaking term A5. 

We would like to mention here that due to the isolated top quark signature, Higgs FCNC 
$ — > tc event may be easy to search for experimentally. However, it is very difficult to 
isolate Higgs FCNC $ — > sb events from the background. 

The paper is organized as follows. In the next section, the 2HDM is introduced. Rel- 
evant couplings are given, theoretical and experimental constraints on 2HDM parameters 
are discussed. In the third section, we will study the effects of 2HDM on Br(Q — > sb) 
which are evaluated in 2HDM-I and 2HDM-II. A comparison between Br(h° — > sb) and 
Br(h° — > 77) is also discussed. Our conclusion is given in section 4. 

2 The 2HDM 

Two Higgs Doublet Models (2HDM) are formed by adding an extra complex SU (2) L <g> 
U(l)y scalar doublet to the SM Lagrangian. Motivations for such a structure include 
CP-violation in the Higgs sector and the fact that some models of dynamical electroweak 
symmetry breaking yield the 2HDM as their low-energy effective theory [17]. 

The most general 2HDM scalar potential which is both SU(2) L <g> U(l)y and CP in- 
variant is given by [18]: 

V(^, $ 2 ) = \ 1 (\$ 1 \ 2 - v\f + A 2 (|$ 2 | 2 - vlf + A 3 ((|$ 1 | 2 - vl) + (|$ 2 | 2 - vl)f + 

A 4 (|$i| 2 |$ 2 | 2 - \^2?) + A 5 (3?($+$ 2 ) - v x v 2 ) 2 + A 6 [3($+$ 2 )] 2 (3) 

where <£>i and <£> 2 have weak hypercharge Y=l, v\ and t> 2 are respectively the vacuum 
expectation values of $1 and $ 2 and the \ are real-valued parameters. Note that this 
potential violates the discrete symmetry $j — > — $j softly by the dimension two term 
\^t{^i^2)- The above scalar potential has 8 independent parameters (Aj)j = i v .. 5 6, i>i and 
f 2 . After electroweak symmetry breaking, the combination v \ + v\ is thus fixed by the 
electroweak scale through v \ + v \ = (2\/2G F )~ 1 . We are left then with 7 independent 
parameters. 

Meanwhile, three of the eight degrees of freedom of the two Higgs doublets correspond to 
the 3 Goldstone bosons (G^, G°) and the remaining five become physical Higgs bosons: H°, 
h° (CP-even), A (CP-odd) and H ± . Their masses are obtained as usual by diagonalizing 
the mass matrix. The presence of charged Higgs bosons will give new contributions to the 
one-loop induced FCNC couplings, as shown in Fig. (1) d u — > d±s. 

It is possible to write the Aj in terms of physical scalar masses, tan j3, a and A5 (see [19] 
for details). We are then free to take as 7 independent parameters (\i)i=i,... t 6 and tan/3 or 
equivalently the four scalar masses, tan/5, a and one of the Aj. In what follows we will 
take A 5 as a free parameter as well as the physical masses and mixing. 

We list hereafter the Feynman rules in the general 2HDM for the trilinear scalar cou- 
plings relevant for our study. They are written in terms of the physical masses, a, f3 and 
the soft breaking term A 5 : 
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H°H + H~ = — — ^— (M^o(cos 3 /3sina + sin 3 (3 cos a) + M%± sin2/3cos(/3 - a) 
My/ sin 2p 

-sin(/3 + a)A 5 t; 2 ) (4) 
H°H+G- = ^L- sin(/3 - a){M^ - M|±) (5) 

h°H + H~ = ~ %9 (M%o (cos a cos 3 (3 - sin a sin 3 /3) + M^± sin 2/3 sin(/3 - a) 

-Miy sin 2/3 

-cos(/3 + a)A 5 v 2 ) (6) 
ft°lf + <T = ^ cos(/3 - a)(M 2 , - Ml ± ) (7) 

A°H+G-= ^f(M| ± -M 2 ) , f= 2 J$L (8) 

We need also the couplings of scalar boson to a pair of fermions both in 2HDM-I and 
2HDM-II. In those couplings, the relevant terms are as follows: 

tftttxM^- , H°ito,M t ^- , A°itoc^- 2HDM - 1, II (9) 
sin p sin p tan p 

h % b(X M b C ^- H o~ bb(xMb s j^ A o bboc J^ 2HDM-I (10) 
sin/3 sin/3 tan/3 

— sin rv — cos rv — 
h°bbocM b , H°bbocM b , A°bb<xM b tan/5 2HDM - II (11) 

cos p cos p 

(H-bt) L oc , (H-bt) R oc 2HDM - I (12) 

v 1 tan/3 v ' tan/3 v y 

(# _ fa) L oc M fe tan/3 , oc 2HDM - II (13) 

tan p 

In this analysis, we take into account the following constraints when the independent 
parameters are varied. From the theoretical point of view: 

i) The contributions to the Sp parameter from the Higgs scalars [20] should not exceed the 
current limits from precision measurements [21]: \Sp\ < 0.001. We stress in passing that 
the extra contribution to Sp constraint [20] vanish when we take M#± = Ma (A4 = Xq). 
Under this constraint the 2HDM scalar potential is 0(4) symmetric [22]. In this case 
(H + , A , H~) form a triplet under the residual global SU{2) of the Higgs potential. It is 
this residual symmetry which ensures that p is equal to unity at tree level. One conclude 
then that large splitting between M H ± and Ma could violate \5p\ < 0.001 constraint. 

ii) From the requirement of perturbativity for the top and bottom Yukawa couplings [23] , 
tan/3 is constrained to lie in the range 0.1 < tan/3 < 70. 

Hi) It has been shown in [24] that for models of the type 2HDM-II, data on b — > 57 imposes 
a lower limit of M H ± > 350 GeV. In type I 2HDM, there is no such a constraint on charged 
Higgs mass [24]. In our numerical analysis we will ignore this constraint in order to localize 
regions in the 2HDM parameters space where the branching ratios are sizeable. 
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iv) Unitarity and perturbativity constraints on scalar parameters: 

It is well known that the unitarity bounds coming from a tree- level analysis [25, 26] put 
severe constraints on all scalar trilinear and quartic couplings. The tree level unitarity 
bounds are derived with the help of the equivalence theorem, which itself is a high-energy 
approximation where it is assumed that the energy scale is much larger than the Z° and 
W ± gauge-boson masses. We will use, instead of unitarity constraints, the perturbativity 
constraints by assuming that all Aj satisfy: 

\K\ < 4vr. (14) 

Those perturbative constraints on the A, allow us to investigate a larger parameter space 
than the one allowed by unitarity constraints. 

We would like to mention also that when performing the scan over the 2HDM parameters 
space, we realize that for some points the widths r$ of the scalar particles become bigger 
than their corresponding masses: T$ > M$ (<3> = h°, H°, A , H^). This happens both 
when we impose tree level unitarity constraints and/or perturbativity constraints. The 
width becomes large specially when the pure scalar decays like H° — > h°h°, H° — > H + H~ , 
h° -> H + H~, H° -> A°A° and h° -> A°A° are open. We find it is natural to add to the 
above constraints the requirement that the width of the scalar particles remains smaller 
than the mass of the corresponding particles: 

T$ < M$ (15) 

From the experimental point of view, the combined null-searches from all four CERN 
LEP collaborations derive the lower limit M H ± > 78.6 GeV (95% CL), a limit which applies 
to all models in which Br(if ± — > tu t )+ Br(iJ ± — > cs)—l. For the neutral Higgs bosons, 
OPAL collaboration has put a limit on h° and A masses of the 2HDM. They conclude 
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Figure 2: SM width and Branching ratio for H — > tc (left) and H — > sb (right) as a function 
of Higgs mass. 



that the regions 1 < M h < 44 GeV and 12 < M A < 56 GeV are excluded at 95% CL 
independent of a and tan/3 [27]. For simplicity we will assume that all scalar particles 
masses are > 90 GeV. 



3 Higgs FCNC in 2HDM 
3.1 Higgs FCNC in SM 

Before presenting our results in 2HDM, we would like to give the Branching ratio of H — > ic 
and H — > sb in the SM. To our best knowledge, the first calculation for Br(H — ► sb) has 
been carried out in [28]. However, in [28], numerical results have been given only for a very 
light Higgs boson Mh = 9 GeV. Recently a new estimation, using dimensional analysis 
and power counting, has appeared both for Br(H —* sb) [14] and Br(H — > ic) [11]. We 
refer the reader to [11, 14] for more details on those estimations. Here we present exact 
result based on diagrammatic calculations both for Br(H — > sb) and Br(H — > ic). We 
give numerical results for the width as well as for the branching ratio. 
The Feynman diagrams contributing to those process in SM are depicted in Fig .(1) di — ► 
dio- In the case of H — > ic, in Fig. (1) (fx, / 2 ) = (t, c) and // = d, s, b, while for H — > sb 
(/i>/2) is (b,s) and j[ = u,c,t. The full loop calculation presented here is done with the 
help of FormCalc [29]. FF and LoopTools packages [30] are used in numerical analysis. 
The numerical results shown in eqs. (1,2) is derived by FormCalc [29]. 

In the SM, as expected, the branching ratio of H — > tc and H — > sb are very suppressed 
due to GIM mechanism. The branching ratio is very small in both cases for higher Higgs 
mass M H > 2M Z where H -> W + W~ and H -> are open. 

Both in SM and 2HDM, the decay widths T| M and T| HDM of scalar particles: $ = H SM , 
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h°, H°, A and ^ are computed at tree level as follows: 

r| M = £r($ - ff) + r($ ^vv) 

f 

r |HDM = r($ -> //) + r($ -> yy) + r($ -> + r($ -> #^) (16) 
/ 

QCD corrections to $ — > // and $ — > {gg, 77, jZ, V*V* , VV* , V*Hi} decays are not 
included in the widths. The decay widths of the Higgs bosons are taken from [31]. 

For a Higgs mass heavier than 250 GeV, we get branching ratio of the order 10~ 14 — > 
KT 12 (resp KT 10 -> 10- 9 ) for H -> tc (resp H sb). 

In the case of H ^ sb, the branching ratio is enhanced for Higgs boson mass of the order 
Mn ~ 100 — > 120 GeV where the width of the Higgs is very narrow. We have plotted in 
Fig. (2) both the decay width and the branching ratios of H — > tc (left plot) and H — > sb 
(right plot) as well as the branching ratio of H — > As it can be seen from the right 

plot Br(H — > sfe) is two orders of magnitude smaller than Br(H — > fi + fi~). 
Since the decay width of — > tc is very suppressed, the threshold for tt production is 
absent in Fig. 2 (left). The situation is slightly different for H — > s6 where the decay width 
of H ^ sb is about 6 order of magnitude larger than decay width of — > tc. From the 
right plot of Fig. (2) one can see that the Br of H — > s6 is smaller once the ti threshold 
has been passed. 

3.2 h° -> s6 

Turning now to the 2HDM Higgs bosons FCNC couplings $ ^ s5, $ = h°,H°,A°. The 
Feynman diagrams are depicted in Fig. (1). The amplitude is sensitive to the <&H + H~ and 
$/J ± G =F couplings through diagrams di2, 13,14 as well as to the $tt and (H~M)l,r couplings 
through diagrams ^11,12,13,14- In 2HDM, it is expected that the dominant contribution to 
the amplitude of $° — > sb comes from diagram dyi- The amplitude of d\i is proportional 
to the trilinear Higgs coupling §°H + H~ and is given by ($ = h°, H°): 

M dl2 = VH+H-^ M l Jt\ [(1 + tan/^Q, + d + C 2 ]t)(M s )l±^ M (M b ) (17) 
o7r tan p oivi-^s-^ z 

where we have neglected the strange quark mass. In the conventions of [29] , the arguments 
of the Passarino-Veltman functions Q are {M^, M 2 , M^ ± , M 2 , M^ ± }. The Yukawa 
coupling Yb of the bottom is model dependant and is given by Yj, = — 1/ tan/3 (resp Y& = 
tan/5) for 2HDM-I (resp 2HDM-II). 

In 2HDM-I, 1 + tan/3F b = 0, the amplitude of d 12 is enhanced by tan |^ factor for small 
tan/5 as well as by the trilinear coupling §°H + H~. 

The diagram d n is sensitive to the coupling $°tt. It is clear from equation (9) that the 
top effect is enhanced for small tan f3 in the case of CP-odd A boson. While in the case of 
CP-even H° and h°, the couplings are enhanced at small tan/3 and large sin a (resp large 
cos a) for H° (resp h°). Consequently, our numerics are presented for small tan (3 = 0.3, 
sin a = 0.1 for h° — > sb and sin a = 0.95 for H° — > sb. 
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Figure 3: Contours for 2 x Br(h° -> sb) in 2HDM-II tan/5 = 0.3 (left), tan/3 = 1.5 (right) 
in the (M H ±, A 5 ) plane with M h = 110 GeV, M H = 180 GeV, sin a = 0.1 and M A o = M H ± 

We also give other numerical results for specific 2HDM parameters where Br(h° — > sb) 
and Br(H° — > sb) get their maximum values without violating 5p and perturbativity 
constraints. 

We show in Fig. (3) contour plots for Br(h° -> sb) in 2HDM-II tan/3 = 0.3 (left) 
and tan/3 = 1.5 GeV (right) in the (Mh±, A5) plane. A5 is varied in the perturbative 
range \\$\ < 4ir. The other inputs are = 110 GeV, Mh = 180 GeV, sin a = 0.1 and 
Mao = M H ±. The width T h o is computed at tree level according to eq. (16). Since the 
mass of h° is taken at 110 GeV, only light fermions contribute to T^o and so the width 
is very narrow and is of the order 57 x 10~ 4 (resp 83 x 10 -5 GeV) at tan/3 = 0.3 (resp 
tan/3 = 1.5). Such narrow width could enhance the branching ratio Br(h° — > sb). We 
would like to mention first that for this set of parameters, the perturbativity of scalar 
quartic couplings Aj is violated around A5 > 5.5. We get |Ai| > 47r for tan/3 = 0.3, while 
for tan/3 = 1.5 there is no such bound. 

Large branching ratios can be obtained for light charged Higgs mass. This can be seen 
in the left panel black and blue areas of Fig. (3) which correspond to small tan (3 = 0.3 and 
large | A5 1 . In those areas the coupling h°H + H~ gets its largest value (see also Fig. (4)). 
In this case one can obtain branching ratio in the range: 10 -4 < Br(h° — > sb) < 6 x 10~ 4 
for Mh± < 200 GeV , A5 < —1.2 and A5 > 3. For charged Higgs mass greater than 200 
GeV, there is also a region where the branching ratio can be in the range 10~ 5 — > 10~ 4 . 
This can be achieved by taking large and negative A5 < — 1. In the case of positive A5 and 
Mh± 250 GeV, the branching ratio decreases to a value < 10~ 5 . 

When tan/3 = 1.5, the coupling h°tt is reduced, and we are left only with a small 
region where the branching ratio Br(h° —>■ sb) is of the order 10~ 5 — > 10~ 4 for Mh± ^ 250 
GeV and large |A 5 | > 5. In both plots (left and right), the coupling h°H + H~ reaches its 
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Table 1: Maximum Branching ratios of h° — > sb in 2HDM-I and II and corresponding 
2HDM parameters, all masses and decay width are in GeV. In Br and widths T columns, 
the upper row is for 2HDM-I and the down row is for 2HDM-II 



minimal value in the region where A5 ~ — > 2, which explains why the branching ratio is 
so small in this region. 

Now we turn to the case where M H ± ^ M A , 5p 7^ 0. We have performed a system- 
atic scan over the full 2HDM parameters space taking into account 5p and perturbativity 
constraints. The maximum branching ratios found for h° — > sb in 2HDM-I and II are 
displayed in table 1. We show not only width and Br of h° — > sb but also the width and 
Br of h° — > 77 for comparison. The total width of the Higgs r^o is also given. When T h o 
becomes comparable to the width of h° — > sb and/or h° — > 77, those decays widths have 
to be included in the total width I\o in order to compute the Br 5 b and Br 77 . 
The first three columns of table 1 are for 2HDM parameters. From 4th to 8th columns we 
give Br and widths. In those columns, the upper row is for 2HDM-I and the down row is 
for 2HDM-II. 

In 2HDM-I, Br(h° -> sb) of the order IO" 3 can be reached in the limit sin a -> -0.98 
(a — > — 7r/2) and small tan (3 < 0.5. In fact, this limit (a — > —n/2) is very close to 
fermiophobic scenario a = ±7r/2. In the fermiophobic limit, all couplings of h° to down 
quarks and leptons are suppressed eq. (10). In this limit, h°tt is also suppressed eq. (9). The 
width of light Higgs h° {Mh < 160 GeV) is then very tiny in the limit sin a — > —0.98. This 
tiny width together with large h°H + H~ are the sources of enhancement of the Br(/i° — > sb) 
to IO- 3 level. This can be seen in the first, second and third lines of table 1 
In 2HDM-II, the couplings of h° to down quarks and leptons are suppressed for sin a pa 0.1 
eq. (11). Hence, the width of light Higgs (M^ < 160 GeV) is very tiny in the limit 
sin a PS 0.1. Moreover, in this limit, the coupling h°tt is enhanced in both models 2HDM-I 
and II. The decay width T(h° — > sb) which was ~ IO -6 for sin a = —0.98 is of the order 
~ IO -5 for sin a = 0.1. Consequently, the Br(h° — > sb) reaches IO -3 . 
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Figure 4: 2 x Br(h° -> 77) and 2 x Br(h° -> s6) in 2HDM-I and II, m H ± = 100 GeV, 



tan/? = 0.3 (left) and tan/? 
Fig. (3). 



5 (right). All the other parameters are the same as in 



In this scenario, as one can see from table 1, the Br(h° — > 77) in 2HDM-I is 2 x 10 -4 which 
is smaller than Br(/i° — > sb) = 9 x 10~ 4 . This is mainly due to the fact that the trilinear 
coupling h°H + H~ is very suppressed in this scenario (see more details in next section). 
As one can see from the last line of the table 1, there exist also values of sin a = 0.18, far 
from fermiophobic scenario but with small A5 = 0, where Br(/i° — ► sb) can be of the order 
10~ 3 . 



3.3 Can h° — > sb compete with h° — » 77? 

It is well known that the decay h° — > 77 is loop induced and so is suppressed. In the 
SM, the branching ratio Br(H SM — > 77) is about ~ 10~ 3 for Higgs mass in the range 
Mjj = 100 — > 160 GeV. Hence, with maximum branching ratio for h° — > sb of the order 
1 x 10~ 4 6 x 10" 4 in 2HDM-I or II, it is legitimate to compare h° — > 77 and h° — > sb in 
2HDM-I or II. Of course, even if h° — > sb and h° — > 77 has a competitive branching ratio, 
we should keep in mind that h° — > 77 decay has a clear signature while the FCNC decay 
h° — > sb has not. 

We illustrate in Fig. (4) the branching ratio for h — > sb and h — > 77 both in 2HDM-I 
and II. The charged Higgs mass is fixed to 100 GeV. It is clear that in the case tan/3 = 0.3 
h° — > 77 is about one order of magnitude bigger than h° — > sb. While, in the case 
of tan/? = 5 h° — > 77 is more than four orders of magnitude bigger than h° — > s6. 
This is because at tan/3 = 0.3 (resp tan/? = 5) the W loop are suppressed by a factor 
h W + W~ oc sin(/? - a) w 0.2 (resp enhanced by /i W + W - oc sin(/? - a) w 0.96). All the 
dips observed in the plots correspond to the minimum of the coupling h°H~H + . Those dips 
are not located at the same A5, this is due to a destructive interference with others diagrams. 
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Figure 5: Contours for 2 x Br(H° -> sb) in 2HDM-II in the plan (M H ±, A 5 ) M H = 140 
GeV (left), (M fl±1 M#) A 5 = 5 (right) with tan/5 = 0.3, M h = 110 GeV, sin a = 0.95 and 
M A o = M H ± 



When h°H H + coupling is very suppressed, it may be possible that the Br(/i° — > sb) could 
be higher than Br(h° — > 77) as it can be seen both in the left plot of Fig. (4) for A5 = 2.5 
and in table 1 for sin a = 0.1 in 2HDM-II. 

However, even if Br(/i° — ► 77) and Br(/i° — > sb) become comparable, we should keep in 
mind that h° — ► 77 has a very clear signature while h° — > sb does not. 

An interesting feature of the 2HDM-I, is its fermiophobic scenario. The light CP- 
even Higgs h° of the 2HDM-I is fermiophobic in the limit a — > 7r/2, all /i° couplings to 
fermions vanishes for a = tt/2 [32, 18]. If h°, with a mass in the range 100 — > 160 GeV, 
is fermiophobic the dominant decay mode is /i° — > 77. It has been shown in Ref. [33] that 
in the fermiophobic limit, the branching ratio of the one loop induced decay h° — > bb* is 
below 10% — > 30%. As the decay h° — > sb is concerned, we have checked by systematic 
scan that in the fermiophobic limit, the decay width of h° —> 77 is more than one order of 
magnitude bigger than the width of h° —>■ sb. 

3.4 H° sb 

We now discuss the heavy CP-even decay H° — > s6. Our numerical results are shown in 
Fig. (5). To maximize the coupling H°tt, we choose of course small tan/3 w 0.3 and large 
sin a w 0.95. In the right plot of Fig. (5), we show contour plots for Br(H° sb) in the 
plane (M H ±,\ 5 ) for M H = 140 GeV. For CP-even Higgs mass 140 GeV, H° -> W + W~, 
H° -> ZZ, F° tt, F° and ^ifj are not yet open, and so the width is 

narrow. In particular, for the set of parameters fixed here: Mj, = 110 GeV, sin a = 0.95 

*In fact, in the 2HDM, not only the coupling h°"fy and hP^Z [34] can have non decoupling effects, but 
also one loop contribution to h°sb [35] and h°h°h° [36]. 
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Table 2: Maximum Branching ratios of H° — > sb in 2HDM-I and II and corresponding 
2HDM parameters, all masses and decay width are in GeV. In Br and widths T columns, 
the upper row is for 2HDM-I and the down row is for 2HDM-II 



and M A o = M H ±, the width is 73 x 10~ 4 GeV. 

The behavior is similar to what we obtain for Br(h° — > sb). In the black regions (large A 5 ), 
the coupling H°H + H~ is maximal while for As G [0, 2] H°H + H~ is minimal. In the black 
region the branching ratio of Br(H° — > sb) can reach m 7 x 10~ 4 . From the left panel of 
Fig. (5), it is evident that there is a relatively large region in the plane (Mff±,As) where 
the Br(H° -> sb) > 10~ 5 . 

In the right panel of Fig. (5), we show Br(H° — > sb) in the plan (M H) M H ±) for A 5 = 5. 
One can see that when CP-even mass < 2Myy, the decay H° — > W + W~ is not yet 
open. The width T H o is narrow, and so the branching ratio is large. For Mh± < 250 GeV 
and Mh < 2Mw, one can have Br(H° — > sb) > 10~ 5 . Once the CP-even Higgs mass 
Mh > 2M W , the decay H° — > W + W~ is open, and the width is larger than 5 x 10 -2 GeV. 
The Branching ratio Br(H° — > sb) is then reduced. As it can be seen from the right plot, 
for Mh ?t 220 GeV, the branching ratio Br(H° — > sb) is less than < 10 -6 . 

In case of 2HDM-I, both h°tt, H°tt and {H~bt)n couplings are the same as in 2HDM- 
II, while (H~M)l which is proportional to Mj,tan/3 in 2HDM-II is now proportional to 
Mft/tan/3. For small tan/5 ~ 0.3, both Br(h° — > sb) and Br(H° — > sb) are of the same 
order as in 2HDM-II, while for large tan (3 those Branching ratios are less than about 
w 10~ 6 . 

In case where Mh± ^ Ma (Sp ^ 0), we present our results of maximum branching ratios 
of H° -> sb in the table 2. It turns out that in 2HDM-I (resp 2HDM-II), Br(#° -> sb) 
reach 10 -3 for small sin a ~ 0.1 (resp large | sino;| m 0.9). The interpretation is the same 
as in the case of light CP even Higgs h°. In 2HDM-I (resp 2HDM-II), the couplings of H° 
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Figure 6: 2 x Br(A° — > sb) (left) and CP-odd A width T^o (right) as function of Ma for 
several values of tan /3 



to down quarks and leptons are suppressed for sin a ~ 0.1 (resp |sina| ~ 0.9). In those 
cases the total Higgs width is very tiny and so the branching ratio of H° — > sb is enhanced. 
Of course, Br(H° — > sb) reach 10 -3 only for light charged Higgs, which is strongly disfavored 
by b -> S7 constraint [24] in 2HDM-II. 

From table 2, one can see also that in 2HDM-II and for sin a rs {0.9, —0.98} the Br(H° — > 
sb) and Bi(H° — > 77) are of comparable size. This is again mainly due to the suppression 
of the coupling H°H + H~ in those limits. 

As in the case of light CP-even Higgs h°, there exist values of sin a = —0.58 far from 
fermiophobic limit with small A5 = where Br(_ff° — ► sb) can reach 10 -3 . 

3.5 A -> sb 

Let us now look at 2HDM contribution to A — ■> sb. Since A is CP-odd, it does not couple 
to a pair of charged Higgs. The only pure trilinear scalar coupling which contributes 
to A -> sb is A°H ± G^ eq. (8). Unlike the couplings H°H + H~ and h°H + H~ eqs (4,6), 
which depend both on Higgs masses, tan f3 as well as the soft breaking term A5, the coupling 
A°H ± G T depends only on the splitting Mfj ± — M\. As mentioned above, such splitting 
should not be too large, otherwise the 5p constraint is not satisfied. As one can read 
from eqs. (9,13), the couplings A°it and (H~bt)n are proportional to Mt/tan/3. Hence 
enhancement is expected at small tan/?. 

As we stressed before, our 2HDM parameters in this case are: tan/3, Ma and Mh±. 
For simplification, we use the MSSM sum-rules to fix charged Higgs mass and a by using 
tan p, CP-odd mass Ma and a SUSY scale which we take at 1 TeV. CP-odd mass will be 
varied from 100 GeV to 600 GeV without worrying about perturbativity. tan (3 is taken to 
be > 0.1. 

We present our numerical results for A — > sb in 2HDM-II in Fig. (6). As can be seen 
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from the left plot, the Branching ratio Br(A° — > sb) is greater than 10 5 only for small 
tan (3 « 0.1 -> 0.35 and light M A and M H ±. For light M A < 200 GeV and low tan (5 < 1, 
the width of A is still small and so the branching ratio is enhanced. For Ma ?t 200 GeV, 
the decay A — > is open and the decay width T^o increases. Therefore, the branching 
ratio is reduced. Note that for tan/5 = 0.2, we cut off the curve at Ma ~ 400 GeV where 
the width T^o starts to be greater than M A - At large tan/5, due to the bottom Yukawa 
coupling, both the partial width T(A° — > sb) and total width T^o are enhanced, and the 
branching ratio is saturated in the range [10~ 6 , 10 -5 ]. 
The situation is almost the same in 2HDM-I. 

4 Conclusions 

In the framework of the 2HDM with natural flavor conservation, we have studied various 
Higgs FCNC $ — > sb. The study has been carried out taking into account the experimental 
constraint on the p parameter and also perturbativity constraints on all the scalar quartic 
couplings Aj. Numerical results for the branching ratios have been discussed. We empha- 
sized the effect coming from both top and bottom Yukawa couplings and pure trilinear 
scalar couplings such as h°H + H~~ and H°H + H~. 

We have shown that, in 2HDM-I and 2HDM-II, the branching ratios of Higgs FCNC 
{h°, H°, A } — > sb are enhanced to the range of 10~ 4 — >• 7x 10" 4 for small tan /5, rather light 
charged Higgs boson and large soft breaking term A5. The branching ratio of Br({h°, H } — > 
sb) can be pushed to 10~ 3 level when sin a is close to fermiophobic limit (sin a; ~ —0.98) 
or sin a: m 0.1 and even for sin a far from those limits but with small A5 = 0. 
Charged Higgs mass of 2HDM-I is not constrained by b — > S7, Br({h°, H } — > sb) can be 
of the order 10~ 4 — > 10~ 3 for light charged Higgs which is comparable to size of SUSY pre- 
dictions [12, 14]. Those branching ratios rates, could still leads to large number of events 
at LHC [11]. 

In 2HDM-II with b — > 57 constraint, branching ratios of {h°, H } — > sb are smaller than 
10~ 5 (resp 10~ 4 ) for tan/5 > 1 (resp tan/5 < 1). 

In the case of light CP-even m h o ks 100 — >• 160 GeV, we have also shown that the branching 
ratio of Br(h° — > sb) is well below Br(h° — * 77) in most of the case. This is also the case 
in the fermiophobic scenario of 2HDM-I. 

One interesting scenario is that both Br(h° — > 77) and Br(h° — > sb) develop a dips for 
some A 5 (see Fig. 4). Those dips are not located at the same A 5 due to the presence of 
diagrams which contribute to h° — > sb but not to h° — > 77. The dip for Br(h° — > sb) is 
located for A5 = 1 while for Br(h° — > 77) it is located for A5 ~ 2.5. For A5 ~ 2.5, we are 
already away from Br(h° — > sb) dip, the Br(h° — > sb) is slightly higher than Br(h° — > 77). 
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